Development of the spectrum of gamma-ray burst pulses 
influenced by the intrinsic spectral evolution and the curvature 

effect 
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ABSTRACT 

Spectral evolution of gamma-ray burst pulses assumed to arise from emission 
of fireballs is explored. It is found that, due to the curvature effect, the integrated 
flux and peak energy are well related by a power law in the decaying phase of 
pulses, where the index is about 3, being free of the intrinsic emission and the 
Lorentz factor. The spectrum of a pulse in its decaying phase differs slightly 
for different intrinsic spectral evolution patterns, indicating that it is dominated 
by the curvature effect. In the rising phase, the integrated flux keeps increasing 
whilst the peak energy remains unchanged when the intrinsic emission bears an 
unchanged spectrum. Within this phase, the flux decreases with the increasing of 
the peak energy for a hard-to-soft intrinsic spectrum, and for a soft-to-hard-to- 
soft intrinsic spectrum, the flux generally increases with the increasing of the peak 
energy. An intrinsic soft-to-hard-to-soft spectral evolution within a co-moving 
pulse would give rise to a pulse-like evolutionary curve for the peak energy. 

Subject headings: gamma-ray bursts, gamma-rays, relativity 

1. Introduction 

Recent Swift observations of early X-ray afterglows of gamma-ray bursts (GRBs) re- 
vealed that the previous interpretation of the GRB afterglow, which is based on the mecha- 
nism of external shocks, is unlikely to be the mechanism accounting for these new findings. 
Instead, the internal shock origin and the so-called "curvature effect" might be responsible 
for many of the early X-ray afterglows [1-8]. The curvature effect is a combined effect that 
includes the delay of time and the shifting of the intrinsic spectrum as well as other relevant 
factors of an expanding fireball (see [9] for a detailed explanation). The effect was widely 
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investigated in recent years in both the prompt emission and early X-ray afterglow of GRBs 



Early study of the spectral evolution of GRBs revealed that the general evolution of the 
spectra of these objects is hard-to-soft (e.g., [24]). This issue was also frequently investigated 
in recent years [11], [25-28]. Reference [14] showed that the hardness-intensity correlation 
could be accounted for by the curvature effect. Authors of [9] explored the evolutionary curve 
of the hardness ratio (HRC) of GRB pulses and concluded that the main characteristics of 
the HRC of some GRB FRED pulses are in agreement with what predicted by the curvature 
effect. In addition, they showed that the curvature effect alone could not completely explain 
the observed HRCs. Instead, a hard-to-soft intrinsic spectral evolution might exist and 
probably play a role in producing the observed HRCs. 

Here, we investigate how the intrinsic spectral evolution and the curvature effect combine 
in producing an observed spectrum arising from a fireball pulse. The paper is organized as 
follows. In Section 2, we present formulas that can directly applied to the case of the 
intrinsic spectral evolution. Spectral evolution within a pulse of fireballs radiated with a 
typical intrinsic emission form is investigated in various aspects in Section 3. In Section 4, 
other intrinsic emission forms are considered. In Section 5, we pay our attention to various 
situations of the soft-to-hard-to-soft intrinsic emission. Influence of the Lorentz factor on the 
spectral evolution is investigated in Section 6. Presented in Section 7 is a signature of the 
curvature effect shown by the evolutionary curve of the peak energy in the decaying phase 
of pulses. Conclusions are given in Section 8. 

2. General formulas of flux and count rates for expanding fireballs 

Here we present more general formulas compared with those derived in [13] and [19] so 
that they could be directly applied to various cases of the intrinsic spectral evolution pattern. 
For more general applications, formulas with various variables are provided. 



For a radiation associated with an intrinsic emission independent of direction, the ex- 
pected flux of a fireball expanding with a constant Lorentz factor is [13] 



[4], [7-9], [10-23]. 



2.1. Formulas in terms of the integral of 6 




(1) 



- 3- 



where I v (t', v, vq) is the observer frame intensity; v is the observation frequency; u is the rest 
frame frequency which emits from differential surface <Isq and gives rise to v (see equation 
9 presented below for the relation between I v and z/ ); t is the observation time; D is the 
distance of the fireball to the observer; 9 is the angle, of dsg, of the fireball, to the line of 
sight; t' is the emission time (in the observer frame), called local time, of photons which 
emit from dse] R(t') is the radius of the fireball measured at time t'. The integral range of 
9, 9 min and 6* max , will be determined by the concerned area of the fireball surface as well as 
the emission ranges of the frequency and the local time. As shown in [13] and [19], t' and t 
are related by 

, t - t c - - + cos e 

t' = \ C R \ + t c , 2 

1 — p cos 9 

where t c and R c are constants. Assume that the area of the fireball surface concerned is 
confined within 

9mm 5: 9 < $max (3) 

and the emission time t' is confined within 

tc ^ t m in — t — ^max) (4) 

and besides them there are no other constraints to the integral limit of (1). According to (3) 
and (4), one can verify that the lower and upper integral limits of (1) could be determined 
by 

9 min = cos 1 min{cos6» min , — ma ^ c R J (5) 

and 

t' ■ - t + - 

9 max = cos 1 max{cos6» max , — min c }, (6) 

(Cin - QP + f 

respectively (for a detailed derivation, one could refer to [13] and [19]. Applying the relation 
between the radius of the fireball and the observation time we come to the following form of 
flux (see [13] and [19]): 

2n[(t-t c -^)(3c + R c ] 2 [°™ I v (t,v,v )cos8sm8 JQ 



From (2) one finds that, for any given values of t and 9, local time t' would be uniquely 
determined. If this value of t' is within the range of (4), then there will be photons emitted 
at i! from the small surface area of 9 reaching the observer at t [when 9 is within the range 
of (3), this small area would be included in the above integral, otherwise it would not]. 



-4- 



Obviously, for a certain value of 9, the range of t depends on the range of t'. Inserting (2) 
into (4) and applying (3) we obtain 



(1-/3 cos min )Cin + (tcP ~ ^) cos 9 min + f < t 

< (1 - COS max )Cax + - f ) cos 9 max + f U 



It suggests clearly that observation time t is limited when emission time t' is limited. 

Let to be the proper time corresponding to the coordinate time t' and z/ be the co- 
moving frequency corresponding to v. The latter two are related by the Doppler effect. It is 
well known that the observer frame intensity I u {t', u, i>o) is related to the rest frame intensity 
loAtoi u o) by 

/„(*>, i/o) = (-) 3 / ,„(to,^o), (9) 

which could be written by 

J - (t ^^ o)= (i-/3cos^r3 (10) 

when applying the Doppler effect. The flux is now able to be written by 

2,r[(i-i c --S)/5c+BJ 2 [»"•* /„,„(«„, v ) cos sin 9 
/.(«) = ^ ^ (1 _ ;3cose)5 <». (") 

where # m j n and ^ m£H are determined by (5) and (6), respectively, z/o and v are related by the 
Doppler effect, and t is confined by (8). 

As to and i! represent the same moment, according to the theory of special relativity, 
they are related by t' — t c — T{t — t ,c), where t 0iC is a constant (here we assign t' = t c when 
to = t , c ). Assign i! = t' m J t()=to m . n and t' = t max | t()= i 0max . We get t min = r(t , mi „ - t 0)C ) + t c , 
and t max = r(t 0imax — t ,c) + t c - Then the constraint of (4) is identical to 

^0,c < ^0,min < ^0 ^ ^O.max- (12) 

We then get from (5)- (6) that 



0mm = cos 1 min{cos6> min , 0,m ™ °' c - ^ o Rc — — } (13) 



(^0,max ~ ^0,c)r + ^ — (t — t c ) 
(^0,max ~ t 0iC )rP + ^ 

and 

n -1 r n (*0,min — to,c)F + — — (t — t c ) 

#max = cos max{cos# max , c p }. (14) 



In addition, from (2) and t' — t c — T(t — to, c ) we obtain 



_ t - t c - £ + ^£ cos 

to= (l-^cos^)r +to - c (15) 
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and from t' mhl = r(£ ,min - *o,c) + t c , t' max = r(t ,max - t , c ) + t c , and (8) we gain 

(1 - /3cOS0 min )[(* o ,min - t 0tC )T + t c ] + (t c f3 - COS min + f < t 



< (1 - /3cOS# max )[(t ,max ~ t 0) c)r + * c ] + (*c/3 - ^) COS # max + 



(16) 



c 



When calculating the flux from (11), we need to convert variable t to 9 via (15) and 
convert u to 9 according to the Doppler effect u — (1 — /3cos#)rz/. 

Light curves of gamma-ray bursts are always presented in terms of count rates within 
an energy range. The count rate within energy channel [i/i, vq\ is determined by 

dn(t) _ r Mt\ dv (17) 



dt J Vl hv 



Inserting (11) leads to 

dn(t) 2n[(t-t c - ^)(3c + R c ] 2 [*"•** p I 0>u (t ,u ) , . costfsin 



dt - in*h h_ X, - %- u^o/ (K (18) 

In the case that 9 could be treated as a constant, from u — (1 — /3cos6>)IV one gets 
^| e)=const = ^a. Replacing variable z/ with u we get from (18) that 

dn(t) _ 2n[(t-t c - ^)[3c + R c ] 2 p- 2 I 0jl/ (t Q , v ) costfsinfl 



n> 2 V(to^o) cos g sing 



where z/ 0j i = (1 — (3 cos9)Ti>i and u 0i2 = (1 — /9cos^)rz/ 2 . 

2.2. Formulas in terms of the integral of cos 6* 

In terms of the integral of cos#, the forms of the above formulas become simpler. Let 

fi = cos9, (20) 



(^O.max — ^0,c)r + ^ — (t — t c ) 
(^0,max — to ;C )rP + ^ 



AW = cos 6> min = minjcos 6> min , ^ : ^ a , Re } (21) 



and 



(^O.min — ^0,c)r + — — (t — t c ) 

x = max|cos(7 max , 

where (13)-(14) are applied. 



AW = cos 6» max = max{cos 6» max , — ■ — : c fl }, (22) 

(f0,min _ H),cJJ- P + 
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The flux now can be calculated by 

2n[(t -t c - f )(3c + R c } 2 r*™ I 0:U (t , u )fi 

A' min 

where t and \i are related by 



_ 27T[(t-t c -f)/jc + ^ c ] 2 Io tV (to,U )fl J 

m ~ " L in * (23) 



t — t — — + — u 

'° = (i-k H< (24) 

and u and are related by vq — (1 — fln)Yv. 
The count rate is now determined by 

dn(f) 2vr[(t - t c - f )/3c + R c } 2 . p- 2 7 0)I/ (t , ^o) , , ^ 



eft D 2 T 3 h 



'Mmin •'"CI 

where i/q,i = (1 — /3//)IVi and i/ ,2 = (1 _ /3/i)rz/ 2 . 



-u, (25) 



2.3. Formulas in terms of the integral of the proper emission time 

Since the rest frame intensity Io tU (t , u ) is always provided in the form of the function 
of the proper emission time to, it would be more convenient to compute the flux or the count 
rate with formulas in terms of the integral of to. 

From (24) we get 
and from z/ = (1 — Pfi)Tu we find 

f - [g - (t - t c )]( 3 
+ (t - to, c )V(3 



UQ = c_ Vc_ V_ (27) 



According to (26) and (27) one can check that ^ > 0. 

Replacing /i with t , we obtain the following form of formula for calculating the flux 

2ttc 2 4° m 7 V(*o, ^[fo ~ *o,c)r + ?-(*- ic)][f + (to - to,c)I73] 2 dto 

= — D 2 r2{ ^_ [? _ (t _ tcM 2 > ( 2§ ) 

where to.min and to,max are determined by 

t - t c ~ § + v cos 6> max 

^0,min = niax{to,min, 7q f 7 rp; 1" ^0,c} (29) 

(i -/5cos^ max )r 
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and 

t - t c - f + ^ cos 6> m in 
t(),max = min{t ,max, 7Z f 7 — 77; 1" t , c }, (30) 

(1 - (3 cos 9 min )T 
respectively, and z/ and t are related by (27). 

The count rate is now written in the form 

dn(t) _ 2 ^ 2 jfriC J °-^ to ' vHfto ~ *o,c)r + g - (t - Q][g + (t - W)T(3fdt 
dt ~ ' D*I*h{% -[f -(t-t c )]/3} 2 

(31) 

According to (27), in the case that t and to could be treated as constants we get ^-\t =const,t=const 
Replacing variable v with u we get from (31) that 



dnjt) _ 2^ 2 tZli:: ^t^^o][(to - t , c )r + g - (t - t c )][f + (to - to, c )r/^t 

dt £> 2 r 2 M^-[?-(*-*cM 2 

(32) 

where z/ ,i = [f -(? -t+t c )/3]r^/[f + (t -t , c )VP] and z/ ,2 = [f - (f -t + t c )/3]rz/ 2 /[f + 

(t - v c )r/3]. 



2.4. Formulas in terms of the integral of relative timescales 

Assign [19] 

f ' — t t' . — t t' —t 

T — R c 1 'min — R c J 'max — R c ) 



c c c 



_ to — to,c _ to,min — to,c to,max — to,c 

T = r^ ) r 0,min = ^ , r 0,max = , {^) 

c c c 

and 

t - t c - - + — 
r ee ^ (35) 



C 



One finds that 



U ° = l + (3Tr ^ (36) 

r- (1 -cosfl) 

T = l-/3cos£ ' (37) 

< r min <r'< r max , (38) 

and 

1 - cos 6 min + (1 - (3 cos 6» min )r min < t < 1 - cos 6> max + (1 - /3 cos max )r4 ax (39) 
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[which is the range of r within which the radiation defined within (3) and (4) is observable] . 
In addition, one gets r 0imin = ^p, T , max = %S and 

1 - cos 9 min + (l-p cos 6> min ) rr 0imi n < r < 1 - cos 6> max + (1 - /3 cos 6> max )rT , max . (40) 
Let 

^(to, ^o) = hA^o), ( 41 ) 

We get from (28) that 

t ^ _ 27rg V(r , */ )(l - r + r r)(l + r/jr ) 2 rfr 

]v[t) ~ ~ D^{l-f3 + f3rf ' (42) 

where z/ is related to r and z/ by (36), 

~ ^0,min — t c T— 1 + COS # max 

r °- in = ^ = max ^-> ( i-^cos^ max )r } ' (43) 



c 



and 



~ _ ^0,max — to,c _ . r T — 1 + COS 6 min 

r , max - m - mmir , max , ^—^—^—^y. 



(44) 



Let 



We then get 



C(r) = (45) 
dr 



2,Rl Jgjgf ^ggg^ " T + rr ° )(1 + ^ To)2dro 
C(T) = hcD^l-P + Pry • (46) 

where z/ 0j i = ^j^0^~ u i an d ^0,2 = u 2- Note that (42) and (46) hold only when 

(40) is satisfied. In the range beyond (40), f v (r) and C{r) will become zero. 

Formula (46) shows that the profile of count rates of a fireball source is a function of r. 
It is independent of the real time scale t, and independent of the real size, R c , of the source. 
In other words, no matter how large is the fireball concerned and how large is the observed 
timescale concerned, for the profile of the count rate, only the ratio of the latter to the time 
scale of the initial radius of the fireball plays a role. 

Applying (33) and (34) to t' — t c — T(t — £ 0>c ) we get r' = Tt . Substituting it into (37) 
comes to 

r-(l-cosfl) 

r ° = (i-/w)r (47) 
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Generally, observation time r is related to the co-moving time To of the surface with 9 by 

T = (I- (3 cos 9)Tt + I - cos 9. (48) 

For the same emission time To, photons emitted from surfaces with different line-of-sight 
angles would reach the observer at different times. In units of the time scale of the initial 
radius R c , the time delay of photons emitted from 9 relative to that of those emitted from 
9 = is 1 — cos 9. The contraction factor of the emission time is (1 — /? cos#)r which differs 
from surface to surface. Taking 9 = we gain 

t = (1-P)Tt (9 = 0). (49) 

This is the relation that relates the observation time of photons which are emitted from 
9 = with the corresponding co-moving emission time. The time is contracted by a factor 

of (i-/3)r. 



3. Spectral evolution within the period of fireball pulses coming from a 

typical intrinsic emission form 

As revealed in [9], the observed hardness ratio is seen to be harder at the beginning of 
some GRB pulses than what the pure curvature effect (when no intrinsic spectral evolution is 
involved) could predict and be softer at later times of the pulses, causing a so-called "harder- 
leading" problem. An economic mechanism was suggested to account for this problem, which 
is the hard-to-soft evolution pattern of the rest frame spectrum of fireballs (see [9]). In a 
more realistic situation, before the hardest spectrum appears, there might be a short period 
within which the co-moving spectrum undergoes a soft-to-hard phase. This would happen 
when the inner shell possesses a relatively harder core and surrounding the core are attached 
with less dense materials. The hardest spectrum will appear when the core strikes the outer 
shell. Before that the spectrum radiated from the outer shell is softer since the shell is hit 
by relatively less dense materials and then gains a much smaller acceleration and a smaller 
speed. 

Let us investigate the spectral evolution in a simple case where the emission of a co- 
moving pulse with an exponential rise and exponential decay and with a flexible Comptonized 
radiation form from the whole fireball surface is observed. For the sake of comparison, we 
consider three co-moving pulses (or rest frame pulses) with their spectra being unchanged, 
hard-to-soft, and soft-to-hard-to-soft respectively during the period concerned. The co- 
moving pulse with an unchanged spectrum is assumed to be 

7 i \ t i+«r / / u exp(- r "'" T " ) (r min < t < T ) 

V(r , u ) = I u + c exp(-z/o/z/ , c ){ ^ T0 _% )0 ^ ' '. (50) 
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The co- moving pulse (50) is intrinsically the same as that adopted in [9]. [Note that, r' = Tr , 
or t oc t' ] The co-moving pulse possessing a hard-to-soft spectral evolution is assumed to 
be 

7 I \ T l+a c \ I I Ml Mr exp(- r "'° T ° ) (r 0l min < T < Tqq) 

V( r o,^o) = V + c exp[-(r /ro,o)(^o/^o,c)]{ , T (1 -\<n < 'n> 

ex P( —) (to,o < r J 

(51) 

where To,o 7^ 0. The co-moving pulse with a soft-to-hard-to-soft spectrum is assumed to be 

J ( Nr l+a c r ex P[-( r 0,0/T-0)K/^0,c)]exp(-^-^) (r , m i„ < T < T , ) 

° M °' ° J ° ° 1 exp[-(ro/r 0l o)(^Ko)] exp(-2=2M) ^ < Tq) , 

(52) 

where To > is required (as shown below, this will easily be satisfied). As adopted previously, 
we assign r = 10oy + r 0)m i n . According to (4), we assign t' > t c , which leads to r 0)m i n > 0. 
As long as this condition is satisfied, r 0)m i n could be freely chosen. [Note that, in the case 
that t' > t c does not hold, r 0)m i n is also constrained: To, m in > — (/^r) -1 , when (5 > (see 
[19]). Without any loss of generality we take r 0) min = 0. To meet the condition that the 
soft-to-hard period is much shorter than the hard-to-soft phase in the third co-moving pulse, 
we take a d = 2a r and this will be applied to all the three co-moving pulses in the following 
analysis. 

Shown in Fig. 1 are the spectral evolutions of the three co-moving pulses (50), (51), 
and (52), where the co-moving times concerned are r (l) = t 0j o — 2<r r , r (2) = r 00 — ay, 
To (3) = 7b,o, r (4) = r 0i0 + <r d , r (5) = r 0j0 + 2a d , r (6) = r 0j0 + 4a d , r (7) = r 0j0 + 6a d , 
To (8) = r , + 8cr d , r (9) = r 0i0 + 10cr rf , and r (10) = r 0i0 + 12cr d , respectively. 

Let us check how these co-moving pulses give rise to the spectra of fireballs. According 
to (49), the observation times of photons which are emitted from the tip of the fireball 
(i.e., 9 = 0) at the above co-moving times are r(l) = (1 — /5)r(ro,o — 2oy), r(2) = (1 — 
/3)r(r 0l0 - <r r ), t(3) = (1 - /3)r(r 0l0 ), r(4) = (1 - ^)r(r , + <T d ), r(5) = (1 - /3)r(r 0)0 + 2<r d ), 
r(6) = (1 - /3)r(r ,o + 4a d ), r(7) = (1 - /3)r(r 0)0 + 6<7 d ), r(8) = (1 - /3)r(r 0)0 + 8<7 d ), 
r(9) = (1 — /3)r(r 0; o + 10<r d ), and r(10) = (1 — /3)r(r 0j0 + 12cr d ), respectively. Positions of 
these times in the light curves observed by distant observers are shown in Fig. 2, where, 
(46) is employed to calculate the corresponding light curves. Note that, equation (21) in [19] 
could only be applied to a simple case when the time component and the spectral component 
of an intrinsic pulse could be separated (see equation (10) in [19]). The advantage of equation 
(46) is that it could be directly applied to a co-moving radiation which can take any forms. 
For the same reason, (42) will be applied to calculate the observed spectrum in the following 
analysis. 



In Fig. 2, two typical widthes of co-moving pulses are adopted to show how the co- 
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Fig. 1. — Spectral evolutions of co-moving pulses (50) (the lower panel), (51) (the mid panel), 
and (52) (the upper panel), where we take a r = 0.01, = 2<r r , a o c - = —0.6, z/ ,c — IkeVh^ 1 , 
and Iq — 1. Solid lines from the bottom to the top stand for the co-moving times To(l), 
ro(2), and r (3), respectively, and dashed lines from the top to the bottom represent r (4), 
To (5), To (6), t (7), r (8), r (9) , and r (10), respectively. 
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Fig. 2. — Positions of the observation times corresponding to various co-moving times for 
photons emitted from the tip of the fireball (9 = 0) marked in the light curves expected by 
distant observers. The light curves associated with co-moving pulses (50) (lower panels), 
(51) (mid panels), and (52) (upper panels) are calculated with equation (46), where we take 
a d = 2a r , a 0iC = -0.6, v ,c = IkeVh' 1 , T = 100, and C = 27iR 3 J /hcD 2 = 1. The 
width in the rising part of the co-moving pulses is taken as u r = 0.0001 (left panels) and 
u r = 0.01 (right panels), which correspond to "narrow" and "broad" pulses, respectively (see 
[17] and [19]). Solid lines stand for the light curves expected in the third channel of BATSE 
(55keV < E < HOkeV). Dashed lines from the left to the right represent observation times 
r(l), r(2), r(3), r(4), r(5), r(6), r(7), r(8), r(9), and r(10), respectively. 
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moving width influences the observed spectrum. They are ay = 0.0001 and oy = 0.01 in co- 
moving pulses (50)- (52), which correspond to relatively narrow and broad pulses respectively 
(see [17] and [19]). The words "narrow" and "broad" refer to the profile of the observed 
pulses. A "narrow" pulse comes from a local pulse with a small ratio of its width to the 
radius of the fireball [or, the ratio of the co-moving pulse width to the product of the radius 
of the fireball and the Lorentz factor is small; see equation (49)]. A character of this kind 
of pulse is that the pulse possesses a more steadily decaying phase. Or more precisely, 
the deviation of its decaying profile to the so-called standard decaying form is mild (for a 
detailed explanation, see [17]). Conversely, a "broad" pulse arises from a local pulse with a 
large ratio of its width to the radius of the fireball. A character of this kind of pulse is that 
its decaying profile possesses a reverse S-feature deviation from the standard decaying curve 
(or the deviation is very obvious) [17]. 



3.1. Evolutionary pattern of the overall spectrum 

Displayed in Fig. 3 are the spectral evolutions within the period of the expected observed 
pulses of fireballs, which are outcomes of the three co-moving pulses (50), (51), and (52). 
We use equation (42) to calculate the spectra measured at the observation times t(1), t(2), 
t(3), t(4), t(5), t(6), t(7), t(8), t(9), and t(10). The resulted overall spectra show a 
general hard-to-soft evolution pattern within the period of the pulses observed. During 
these observation times, the spectral evolution of "narrow" pulses is mild (see Fig. 3 left 
panels), while for "broad" pulses their spectra develop quite rapidly (see Fig. 3 right panels). 
In the decaying part of light curves of "narrow" or "broad" pulses, the spectral evolution 
patterns are quite similar for the three co- moving pulses (see dashed lines in Fig. 3). This 
must be due to the fact that the decaying part of the pulses is dominated by high latitude 
emission of the fireball, where the curvature effect is important. It is then natural that the 
spectral evolutions show a similar trend. In the rising part of the light curves of pulses, the 
situation is different. The spectral evolution pattern of co-moving pulses obviously influences 
the way the observed spectrum evolves. 



3.2. Evolution of the peak energy of the vf v spectrum 

There are two quantities capable of describing the hardness of spectra. One is the so- 
called hardness ratio and the other is the peak energy E p , the energy where the peak of 
uf v is observed. Evolutions of the two quantities could describe in some extent how the 
spectrum evolves. Recently, the hardness ratio curves associated with different situations 
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Fig. 3. — Spectral evolutions of fireball pulses arising from co-moving pulses (50) (lower 
panels), (51) (mid panels), and (52) (upper panels). The spectra are calculated with (42). 
Each panel in this plot corresponds to the same panel in Fig. 2, where the parameters are 
the same. The only difference is that we take F = 2kR 2 c Iq/D 2 = 1. In the left panels, 
solid lines from the bottom to the top stand for the spectra measured at observation times 
t(1), t(2), r(3), t(4), and r(5), respectively, while dashed lines from the top to the bottom 
represent the spectra measured at r(6), t(7), t(8), t(9), and r(10), respectively. In the right 
panels, solid lines from the bottom to the top stand for the spectra measured at t(1), t(2), 
t(3), and t(4), respectively, and dashed lines from the top to the bottom denote the spectra 
measured at r(5), r(6), r(7), r(8), t(9), and r(10), respectively. 
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were discussed [9]. Here, we pay our attention mainly to the evolution of the peak energy 
within the period of fireball pulses which arise from the three co-moving pulses proposed 
above. 

We employ the following integrated flux (see Section 5) to compare the evolution of E p 
with light curves: 



where v\ and v 2 are the lower and upper limits of the energy channel concerned. 

Displayed in Fig. 4 are the developments of the peak energy in various cases considered 
above. We find that the feature of drop-to-rise-to-drop revealed in the hardness ratio curve 
in [9] is observed in the peak energy evolutionary curve and it holds only for the unchanged 
intrinsic spectrum. In the case of intrinsically hard-to-soft spectrum, we observe a continuous 
hard-to-soft pattern, while in the case of intrinsically soft-to-hard-to-soft spectrum, one finds 
a soft-to-hard-to-soft evolution pattern for the observed spectrum. There is a turning point 
in the E p (t) curve (hereafter EPC). After the turning point, the EPCs associated with the 
three co-moving pulses show the same trend of softening in the case of "narrow" pulses. 
Before the turning point, the observed spectra vary significantly in accordant with their 
intrinsic spectral evolution patterns. This suggests that the spectrum observed before the 
turning point must mainly be determined by the intrinsic evolution pattern, while after 
the turning point, it is dominated by the curvature effect. In the "broad" pulse cases, the 
softening of the spectra arising from the un-changed intrinsic spectrum co-moving pulses are 
different from that of the hard-to-soft and soft-to-hard-to-soft intrinsic spectrum co-moving 
pulses, which leads to a noticeable character discussed below. The sharpness of the turning 
feature might become an indicator showing a pulse observed is "narrow" or "broad" . 

3.3. Comparison of the spectrum of the emission of the tip with that arising 



If the fireball moves relativistically towards the observer but not expands, one would 
observed an entirely different evolution pattern of the spectrum within the pulse concerned. 
Compared with this spectrum, is that arising from the expanding fireball harder or softer? 
An investigation on this issue might bring us some interesting information. The spectrum of 
an ejecta moving towards the observer must be the same as that emitted from the tip of an 
expanding fireball surface when the intrinsic spectra and the Lorentz factors are the same. 

The spectrum of the emission from the tip of the fireball surface must merely be a 
blue-shifting of the intrinsic one. Due to the contribution of high latitude emissions, a 




(53) 



from the overall fireball surface 
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Fig. 4. — Evolutionary curves of the peak energy (upper panels) of fireball pulses arising 
from co-moving pulses (50) (dotted lines), (51) (solid lines), and (52) (dashed lines). For 
the sake of comparison, we present in the lower panels the integrated flux F 2+3 over the the 
energy range of the second (55keV < E < HOfceV) and third (HOfceV < E < ?>2QkeV) 
BATSE channels, which are calculated with (42) and (53). Parameters are the same as they 
are in Fig. 2. As adopted in Fig. 2, the width in the rising part of the co-moving pulses 
is taken as ay = 0.0001 (left panels) and ay — 0.01 (right panels), which correspond to 
"narrow" and "broad" pulses, respectively. 
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deviation of the observed spectrum of the whole fireball from that of the tip is expected. If 
the deviation is small within some period of observation time, then one can estimates the 
form of the intrinsic spectrum from the observed one within this time interval, although the 
blue-shifting factor and the amplifying factor (or, the boosting factor) remain unknown. Let 
us investigate how this deviation evolves and what it depends on. 

3.3.1. In the corresponding overall spectra 

The method for comparing the expected tip spectrum with that of the emission from 
the whole fireball surface is straightforward. One can simply shift the co-moving spectra 
shown in Fig. 1 to high energy range where the emission of the tip is expected, and then 
compares them with those presented in Fig. 3. According to the Doppler effect, photons 
emitted from the tip with frequency u would be blue shifted to v — i/ /F(l — (3) when they 
reach the observer. We will use this shifting factor, 1/T(1 — f3), to shift the spectra in Fig. 
1. 

Shown in Figs. 5, 6 and 7 are the comparisons between spectra expected from the 
emission of the tip (which is merely the shifting of the intrinsic spectrum) and the emission 
from the whole fireball surface, made for co-moving pulses (50), (51), and (52), respectively. 
For the sake of comparison, we normalize the shifting spectra in Fig. 1 to the maximums of 
the corresponding spectra in Fig. 3. 

In the case of the unchanged intrinsic spectrum, the observed spectrum in the rising 
part of the light curve for "narrow" pulses is almost the same as that expected from the 
emission of the tip of fireballs (see the left upper panel of Fig. 5), while it deviates slightly 
from the latter for "broad" pulses (see the right upper panel of the figure). In the decaying 
portion, the observed spectra are obviously softer than that expected from the emission of 
the tip of fireballs, and the softening is very significant for "broad" pulses (see the two lower 
panels of Fig. 5). This indicates that if an unchanged spectrum during the rising part of 
pulses is observed, then the corresponding intrinsic spectrum is likely an unchanged one. 

When the evolution of the intrinsic spectrum is hard-to-soft, the observed spectrum is 
hard-to-soft as well. In the rising part of the light curve of "broad" pulses, the spectrum 
is almost the same as that of the emission expected from the tip of fireballs (see the left 
upper panel of Fig. 6), while in the decaying portion of the light curve it deviates slightly 
from the latter (see the left lower panel of the figure). For "narrow" pulses, the spectrum is 
always harder than that of the emission expected from the tip of fireballs, and the deviation 
is obvious and it is very significant in the decaying portion of the light curve. 




log v /keVh" 1 log v Q /keVh 



Fig. 5. — Comparison between spectra expected from the emission of the tip (dashed lines) 
and the emission from the whole fireball surface (solid lines) in cases of "narrow" (left 
panels) and "broad" (right panels) pulses arising from the co-moving pulse with an unchanged 
intrinsic spectrum, that of (50). Solid lines in the left upper and lower panels are the solid 
and dashed lines in the left lower panel of Fig. 3, respectively. Solid lines in the right upper 
and lower panels are the solid and dashed lines in the right lower panel of Fig. 3, respectively. 
Dashed lines denote the shifting lines of the corresponding co-moving spectra presented in 
the lower panel of Fig. 1, which correspond to the co-moving times (see the caption of Fig. 
1) when photons emitted at these times from the tip of the fireball surface reach the observer 
at the times concerned in Fig. 3 (see the captions of Figs. 2 and 3). 
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Fig. 6. — Comparison between spectra expected from the emission of the tip (dashed lines) 
and the emission from the whole fireball surface (solid lines) in cases of "narrow" (left 
panels) and "broad" (right panels) pulses arising from the co-moving pulse with a hard-to- 
soft intrinsic spectrum, that of (51). Solid lines in the left upper and lower panels are the 
solid and dashed lines in the left mid panel of Fig. 3, respectively. Solid lines in the right 
upper and lower panels are the solid and dashed lines in the right mid panel of Fig. 3, 
respectively. Dashed lines denote the shifting lines of the corresponding co-moving spectra 
presented in the mid panel of Fig. 1. (See Fig. 3 for a more detailed explanation.) 
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Fig. 7. — Comparison between spectra expected from the emission of the tip (dashed lines) 
and the emission from the whole fireball surface (solid lines) in cases of "narrow" (left panels) 
and "broad" (right panels) pulses arising from the co-moving pulse with a soft-to-hard-to- 
soft intrinsic spectrum, that of (52). Solid lines in the left upper and lower panels are the 
solid and dashed lines in the left upper panel of Fig. 3, respectively. Solid lines in the right 
upper and lower panels are the solid and dashed lines in the right upper panel of Fig. 3, 
respectively. Dashed lines denote the shifting lines of the corresponding co-moving spectra 
presented in the upper panel of Fig. 1. (See Fig. 3 for a more detailed explanation.) 
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In the case of the soft-to-hard-to-soft intrinsic spectrum, the observed spectra in the 
rising part of the light curve are softer than, but quite close to, that expected from the 
emission of the tip of fireballs. In the decaying phase, they are much harder than the latter 
for "narrow" pulses and are almost the same as the latter for "broad" pulses. 

3.3.2. In terms of the peak energy 

One might notice that the magnitude of the spectrum of the emission expected from 
the tip of fireballs discussed about is not real since emission from other parts of fireballs is 
much larger than that merely from the tip and hence they are not comparable. However, 
the concept of this spectrum is useful since it does represent the real form of emission from 
the tip. After all, it is the form that determines how "hard" is a spectrum. 

Let us investigate how the spectrum arising from the whole fireball surface deviates 
from that of the emission of the tip of fireballs in terms of the peak energy. We define the 
deviation of the peak energy of the spectrum of the whole fireball emission and that of the 
tip as: 

AE P = E p — Ep jt i P , (54) 
where E p>tip is the peak energy of the spectrum of the emission expected from the tip. 

Displayed in Fig. 8 are the relations between AE P and time derived from various 
cases. Conclusions obtained above are reinforced, a) In the case of the unchanged intrinsic 
spectrum, that arising from the whole fireball surface is always softer than the spectrum 
merely coming from the tip (see also Fig. 5). In the decaying phase of the pulses, the 
softening becomes stronger and stronger, b) When the intrinsic spectrum evolves in a simple 
hard-to-soft pattern, the spectrum from the whole fireball surface is always harder than that 
merely coming from the tip (which is beyond our initial expectation) (see also Fig. 6). c) For 
a soft-to-hard-to-soft spectral evolution, the spectrum of the whole fireball surface becomes 
softer than that from the tip during the rising phase of the pulses, and then turns to be 
harder than the latter in the decaying phase (see also Fig. 7). d) For "broad" pulses, the 
hardening in the decaying phase in the latter two cases is very mild and thus the observed 
spectrum could serve as a good estimator of that from the tip emission, e) In the case of 
"narrow" pulses, the hardening in the decaying phase is so strong that one cannot estimate 
the intrinsic spectrum from an observed one. A consequence of this phenomenon is that in 
the decaying phase, "narrow" pulses are generally harder than "broad" ones. This could 
possibly be checked by current observation. 

One might notice that the deviation of the peak energy is sensitive in the rising part 
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Fig. 8. — Evolutions of AE p (upper panels) for "narrow" (left panels) and "broad" (right 
panels) pulses, calculated with (54). Light curves of the integrated flux ^2+3 (lower panels) 
are also presented (see Fig. 4) for a direct comparison. Parameters and symbols are the 
same as they are in Fig. 4. 
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of pulses. In the case of the simple hard-to-soft intrinsic spectral evolution pattern, the 
observed peak energy could become about 300keV larger than that expected from the tip of 
fireballs when the pulses observed are relatively narrow, and it could be about 70keV larger 
than that of the tip when the pulses concerned are relatively broad. In the decaying phase, 
E p could be about 50 — 150keV larger than that of the tip for "narrow" pulses whilst it is 
almost the same as that of the tip for "broad" ones. 

3.4. Spectral evolution revealed in other aspects 

3.4-1- Relation between the flux and peak energy 

An important aspect capable of revealing spectral evolution in GRB pulses is the well- 
known relation between the flux and peak energy. Revealed in Fig. 2 of [29], the observed 
flux F was found to be significantly correlated with the peak energy E p within a single burst 
for a subset of their sample. In their study, F is the integrated flux in a wide energy range 
which was also adopted in [30] . Here we explore the relation between the integrated flux and 
peak energy expected from a fireball pulse. Several cases discussed above are studied on this 
issue. 

The results are shown in Fig. 9, where the adopted energy range confining the integrated 
flux is (1, lOOOO^eV. A linear relation between the logarithm of the two quantities is seen 
in the decaying phase of pulses. In terms of statistics, the relation is called the hardness- 
intensity correlation (HIC) which was noticed by many authors in gamma-ray bursts [24], 
[28-29], [31-32]. There is a turnover feature in the relation between the two quantities. The 
turnover feature shows a hook-like curve which is sensitive to intrinsic spectral evolution 
patterns. An intrinsic soft-to-hard-to-soft spectral evolution corresponds to a semi-linear 
correlation between logF and logE p in the rising phase of pulses, with its slope being larger 
than that in the decaying phase. In the case of an unchanged intrinsic spectrum, the relation 
in the rising phase is a straight line in parallel with the F axis. When the intrinsic spectral 
evolution pattern is a simple hard-to-soft one, the observed integrated flux F would decrease 
with the increasing of E p in the rising phase. 

Relations between the integrated flux and peak energy in the decaying phase of all 
the six pulses discussed above are calculated. They are found to strictly follow a power 
law: F oc E^, with 7 ~ 3. This is in agreement with what was found by [16] in a recent 
investigation. 
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Fig. 9. — Relations between the integrated flux and peak energy in cases of the unchanged 
intrinsic spectrum (dotted lines), hard-to-soft intrinsic spectrum (solid lines), and soft-to- 
hard-to-soft intrinsic spectrum (dashed lines) for "narrow" (the lower panel) and "broad" 
(the upper panel) pulses of fireballs. Thin lines denote the relation in the decaying phase of 
pulses and thick lines stand for that in the rising part of the light curve of pulses. 
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3.4-2. Relevant relations 

The first relevant relation discussed is that between the peak value of the vj v spectrum 
and the peak energy E p , which was previously investigated by [41]. We study all the six 
pulses discussed above and find that the relation is similar to that between the flux and 
peak energy (the figure is very similar to Fig. 9, which is omitted). We find the same power 
law relation {yf v ) v oc in the decaying phase of pulses and approximately the same index 
7 ~ 3. 

The second is the relation between the flux f v at a fixed energy E and the spectral peak 
energy E p . Here we adopt E = 200keV. All the six pulses discussed above are studied. 
With a slight difference, the relation is similar to that between the flux and peak energy (the 
figure is omitted). Although f u increases with the increasing of E p in the decaying phase of 
pulses, the relation is no more a power law. 

The third is the relation between the integrated flux over the energy range adopted 
above (see Fig. 9) and the hardness ratio. The hardness ratio adopted here is defined by 
EHR = F 3 /F 2 , where F 2 and F 3 are the integrated fluxes confined in the energy ranges 
of the second and third BATSE channels, respectively. Shown in Fig. 10 is the relation in 
various cases discussed in Fig. 9. In the decaying phase of the six pulses concerned, the 
relation is not a power law, but still, the integrated flux increases with the increasing of the 
hardness ratio. 



4. Spectral evolution of pulses associated with other intrinsic emission forms 

The real GRB emission might be more complicated than what discussed above. We 
wonder if different intrinsic emission forms could lead to other conclusions. 

Following [16], we adopt an intrinsic broken-power-law (BPL) spectrum in the emission 
of a pulse. This pulse is assumed to possess a linear rise and linear decay phases. In the 
same way we consider three evolutionary patterns which are unchanged, hard-to-soft, and 
soft-to-hard-to-soft respectively. 

The co-moving pulse is assumed to be 

(JfiU+a (v <v ob ) r0 ~ T °- min (r mi „ < T < T ) 

T («■ i, \ — T S Ku »-b J V U — U ' ' \J T- ,0-T , m in V U ' mm U — U ' U/ 

'"' r,, ^ W @ W K*<^ }{ aSS (-0,0 < r < r , m J ■ (55 > 



where u o b is a function of t . Co-moving pulse (81) is a pulse with a linear rise and a linear 
decay, emitting with a BPL radiation form. In the case of an unchanged intrinsic spectrum, 
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Fig. 10. — Relations between the integrated flux and the hardness ratio. The symbols are 
the same as they are in Fig. 9. 
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u 0i b is taken as vqj, = v fi when T 0tmin < t < T^max (where 1/0,0 is a constant). In the case of a 
hard-to-soft intrinsic spectrum, v o b is assumed to decrease linearly with the increasing of r : 
fob — ro '"_ r °' mm Vo when r 0m j n < r < r 0max . In the case of a soft-to-hard-to-soft intrinsic 

' T ^Ojinin ' ' ' 

spectrum, vo,b is assumed to rise linearly and then drop also linearly with the increasing of 

To- »o,b = r ?~ T ^ min ^0,0 when T , m in < r < r , , ^0,6 = ~ ,m **T~ n ^0,0 when r , < r < r , max . 

'0,0 '0,min ' '0,max '0,0 

One might notice that the broken-power-law model could well approximate the Band 
function model [33] in very low and high energy ranges. We accordingly assign the typical 
values of the low and high energy indexes available from the fit of the Band function model 
to BATSE GRB spectra to the BPL model used here: a = -1 and (3 = -2.25 [34-35]. In 
addition, we assign r 0)m i n = and take r 0imax — t 0i0 = 2(r 0i o — T 0jmin ). We take T = 100 so 
that a direct comparison with those results in Section 3 could be made. 

We study the evolutionary curve of the peak energy (EPC) in the case of the broken 
power law emission (55) emitted with unchanged, hard-to-soft, and soft-to-hard-to-soft in- 
trinsic spectra for "narrow" and "broad" pulses of fireballs (the half width in the rising phase 
of the "narrow" co-moving pulse is taken as 0.0001, and that of the "broad" one is 0.01). 
The analysis shows that conclusions drawn from Fig. 4 hold in this situation (the figure 
is omitted). Noticeable differences are: in the case of "narrow" pulses, the EPC arising 
from the hard-to-soft intrinsic spectrum is significantly larger than that arising from the 
unchanged and soft-to-hard-to-soft intrinsic spectra in the decaying phase; in the case of 
"broad" pulses, the EPC arising from the unchanged intrinsic spectrum is larger than that 
arising from the other two intrinsic spectra in the decaying phase and it seems to deviate 
from the common trend observed in other cases. 

Relations between the integrated flux and peak energy for the three intrinsic pulses 
confined by the three Vq^Tq) functions mentioned above are shown in Fig. 11. The turnover 
features shown in Fig. 9 are observed in this figure. The power law relation between the 
two quantities in the decaying phase of pulses holds for all cases concerned here, where the 
index is found to be 7 ~ 3 as well. In addition, we notice that the whole curves of the 
F — E p relation vary significantly with the intrinsic emission, and the ranges of the power 
law relation depends strongly on the emission as well. 

5. Spectral evolution of pulses expected in other situations of the 
soft-to-hard-to-soft intrinsic emission 

In the analysis of the soft-to-hard-to-soft intrinsic emission of co-moving pulses per- 
formed above, the peak energy of the intrinsic emission concerned is that starting from and 
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Fig. 11. — Relations between the integrated flux and peak energy in the case of the broken 
power law emission (55) with an unchanged intrinsic spectrum (dotted lines), hard-to-soft 
intrinsic spectrum (solid lines), and soft-to-hard-to-soft intrinsic spectrum (dashed lines) for 
"narrow" (the lower panel) and "broad" (the upper panel) pulses of fireballs. Parameters 
or their relations are: a = -1, (3 = -2.25, r 0imin = 0, r 0ima x - r ,o = 2(r 0i0 - r , m i n ), and 
T = 100. The half width in the rising phase of the "narrow" co-moving pulse is taken as 
0.0001, and that of the "broad" one is 0.01. The thick and thin lines stand for the same as 
they do in Fig. 9. 



-29- 



ending at very low energy and the hardest intrinsic spectrum occurs at the time when the 
peak flux of the co-moving pulse appears. This might not be true in practice. In fact, a 
co-moving pulse might start to emit at a high energy band and then evolve to the maximum 
and then drop to the minimum which still remains in a high energy band. In particular, 
some pulses might arise from the emission that the hardest intrinsic spectrum appears ahead 
of the peak flux. We are anxious if conclusions obtained above are affected when taking all 
these into account. 

Two co-moving pulses with a skew soft-to-hard-to-soft spectrum, starting to emit at a 
high energy band and ending its emission at a relatively lower band are considered. One is 
a modified Comptonized radiation form (52): 

(Vo.min < To < 7"0,l) 

7 / x_ r i+« C ( ex Ph(W r o)( z/ o/^o,c)]exp(-^-^) (to,i < r < r , p ) 



exp[-(r /r 0iP )(z/o/t / o,c)] ^Pi- 1 ^^) (to* < To < r ,o) 
exp[-(r /r 0> p)(i/o/fo ) c)] exp(-^^) (r 0j0 < r < r 0)2 ) 



(56) 

where tq. p = r 0i o-Ar and Aro represents the offset of the time the hardest spectrum appears 
ahead of that of the peak flux. We assign r 0i i = ro iP /2 and r 0j 2 = 2r 0iP . In applying co-moving 
pulse (56), we take Ar = ay, ay = 0.01 and T = 100, and adopt all other parameters as the 
same of those adopted in Section 3 (see the caption of Fig. 1). 

The other is a revised broken-power-law emission form which is still described by (55), 
with u ofi being determined by: u 0>b = Z |J^z^ff ^0,0 when r , min < r < r , p , is 0jb = 
2 2(ro,ma^-r ~r u o,o when r 0jP < r < r , max (where r 0iP = r 0i0 - Ar ). In applying co-moving 
pulse (55) coupling this relation, we assign r 0jmin = and take r 0)ma x - ^0,0 = 2(r 

,0 7~0,min)) 

a = — 1, /3 = —2.25, and T = 100, as adopted in last section. We consider only a "broad" 
co-moving pulse and therefore take its half width in the rising phase as 0.01. 

For these two co-moving pulses, we explore the peak energy evolutionary curve and 
the relation between the integrated flux and energy. They are shown in Figs. 12 and 13 
respectively. Fig. 12 shows that the skewness leads to a forward shifting of the peak of 
EPC. It makes the offset of the peak of EPC larger when comparing this peak with the time 
position of the peak count or peak flux (see Fig. 4). The EPC in the decaying phase is 
obviously less affected by the factors considered here. It is interesting that the shape of the 
EPC is that of an incomplete pulse. The lack of the very small values of E p in the rising 
phase must be due to the fact that the co-moving pulses are assumed to start to emit at a 
high energy band. This assumption seems to be reasonable, and therefore one could expect 
the EPC of GRB pulses to possess an incomplete pulse shape if these pulses are suffered 
from the curvature effect. Like what shown in the EPC, the power law relation between the 
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Fig. 12. — Evolutionary curves of the peak energy associated with co-moving pulses with 
a skew soft-to-hard-to-soft spectrum. In the upper panel, the solid line is the curve arising 
from co-moving pulse (56); the dotted line is the dashed line in the upper left panel of Fig. 4. 
In the lower panel, the solid line is the curve arising from co-moving pulse (55) coupling the 
z/ o,6( r o) function presented in Section 5; the dotted line is the curve arising from co-moving 
pulse (55) coupling the soft-to-hard-to-soft vo,b( T o) function presented in Section 4. The 
pulses concerned are "broad" ones (see Section 5 for the details of the adopted parameters). 
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Fig. 13. — Relation between the integrated flux and peak energy associated with co-moving 
pulses with a skew soft-to-hard-to-soft spectrum. In the upper panel, the solid line is the 
curve arising from co- moving pulse (56); the dotted lines represent all the lines in the upper 
panel of Fig. 9. In the lower panel, the solid line is the curve arising from co- moving pulse 
(55) coupling with the uq^tq) function presented in Section 5; the dotted lines denote all 
the lines in the upper panel of Fig. 11. Parameters are the same as thy are in Fig. 12. The 
thin solid line stands for the decaying phase of pulses and the thick solid line denotes the 
rising phase. 
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flux and peak energy in the decaying phase is stubborn (see Fig. 13). The factors concerned 
have no effects on this relation. On the contrary, the relation in the rising phase varies 
significantly according to the evolutionary pattern of intrinsic emission. As shown in Fig. 
13, although the curves in the rising phase are so different, they are all on the right-hand-side 
of the decaying curve, which makes the relation within the whole pulse interval a hook-like 
curve. 

As proposed in Section 3, it will be natural if a co-moving pulse experiencing a soft-to- 
hard spectral evolution phase (which might be quite short) and then a hard-to-soft phase 
during a shock. One thus can expect that the features of the curves studied in this section 
might be common in the GRB pulses if they do arise from the emission of a relativistically 
expanding fireball. 

6. Influence of the Lorentz factor 

In the above analysis, one important factor is not taken into account, which is the 
Lorentz factor of the expanding motion of fireballs. Would it give rise to entirely different 
results? We explore this issue by adopting various values of the Lorentz factor and then 
comparing the results. The peak energy evolutionary curve and the relation between the 
integrated flux and peak energy are investigated. 

We adopt the intrinsic pulses discussed in last section to study this issue, where, param- 
eters other than the Lorentz factor are maintained. The Lorentz factor T = 100 is replaced 
with T = 10 and T = 1000 respectively to show if this quantity plays a role in the two 
relationships concerned. The analysis shows that, although the adopted Lorentz factors are 
significantly different, the features of EPC observed in Fig. 12 are maintained (the figure is 
omitted). Shown in Fig. 14 are the relations between the integrated flux and peak energy, 
calculated with the two Lorentz factors. The shape of the curves of the integrated flux vs. 
peak energy is the same as that shown in Fig. 13. The index of the power law relation 
between the integrated flux and peak energy in the decaying phase changes a little. In the 
case of the revised broken power law emission, the index ranges from 7 = 2.84 (for T = 10) 
to 7 = 2.99 (for T = 1000); and in the case of the modified Comptonized radiation, the index 
changes from 7 = 2.95 (for T = 10) to 7 = 3.07 (for T = 1000). It seems that the larger the 
Lorentz factor, the larger value the index. A noticeable feature revealed in Fig. 14 is the 
shift of the peak energy range, which is expectable due to the Doppler shifting. 
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Fig. 14. — Relation between the integrated flux and peak energy arising from fireballs ex- 
panding with r = 10 (the solid line and the dashed dotted line) and T = 1000 (the dashed 
line and the dotted line) respectively. The dashed dotted line and the dashed line stand 
for the curves associated with the revised broken power law emission, co-moving pulse (55) 
coupling with the ^(to) function presented in Section 5; the solid line and the dotted line 
represent that associated with the modified Comptonized radiation, co-moving pulse (56). 
Parameters other than the Lorentz factor are the same as that adopted in Fig. 13. Thin 
lines stand for the decaying phase of pulses and thick lines denote the rising phase. 
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7. Signature of the curvature effect 

It is known that light curves arising from the emission of an intrinsic 5 function pulse 
with a mono-color radiation over the whole fireball surface share the same profile, a marginal 
decaying curve [17]. Light curves of real emission were found to deviate from this curve by a 
reverse-S feature in their decaying phase and this could serve as a signature of the curvature 
effect. Here we consider a co-moving S function pulse radiated at a fixed energy (a mono-color 
radiation), trying to find if there exists a similar curve in terms of the peak energy. 

The co-moving 5 function pulse with a mono-color radiation is taken as 

IoA T Oi u o) = IoS(t - r 0fi )5(is - i/Qfl) (0 < r ,o). (57) 

Not losing the generality, we take r 0j o = 0. When applying equation (40) we take r , m i n = 
r o,max = since r is confined by r , m i n < r < r 0imax and there is emission only at t = 
t 0j o = 0. We consider the emission from the whole fireball surface and then take 6> min = 
and # max = 7r/2. Thus, from (40) we get < r < 1. The relation between the observed 
frequency and the emitted rest frame frequency is formula (36). Applying (57), the formula 
comes to 

This describes the curve of the development of v that arises from the emission of a co-moving 
5 function pulse radiated with a mono-color spectrum. 

One can check that the maximum of v is reached when r = 0. That leads to z/ max = 
(i-/?)r • Equation (58) then could be written as = 1 _ 1 / ^ 8 /3t - Taking u/u max = 1/2 we get 



1-/3 

Combining these relations we get 



Th = —jf-, which denotes the time when the observed frequency is half of the maximum of v. 



(59) 



Z4nax l+T/T h 

It shows that, in terms of r/r h , the curve z//z/ max is independent of the Lorentz factor. 

Let us study various E p (t) curves discussed above in terms of t'/t' h , where t' is the 
observation time set to the moment when the maximum of E p , E Ptmax , is observed and t' h 
denotes the time t' when the observed E p is half of E Pjinax . [Note that t'/t' h = r/T h since 
the coefficient in relating t' and r is canceled; see equation (35).] As shown in Fig. 4, the 
spectral behavior in terms of E p in the decaying phase of fireball pulses is quite stubborn. 
This is due to the fact that the curvature effect dominates the evolution of the observed 
spectrum in this period. We thus examine haw the evolutionary curve of the peak energy 
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Fig. 15. — Evolutionary curve (the dashed line) of the peak energy E p in the decaying 
phase of a "broad" fireball pulse in the case of the broken power law emission (55) with an 
unchanged intrinsic spectrum (see the caption of Fig. 11). The solid line is the curve of (59). 

in this period differs from the curve of (59) which we call the marginal decaying form of the 
peak energy evolutionary curve. 

We compare the E p {t) curves in the decaying phase of the six pulses discussed in Section 
3 (see the caption of Fig. 2) with the curve of (59) and find that the former well follow the 
latter, where only very minute deviations are observed (the figure is omitted). It is interesting 
that when a deviation is visible, the E p {t) evolutionary curve must possess a reverse-S feature 
relative to the curve of (59). The result indicates that the marginal decaying form (59) could 
serve as a signature of the curvature effect. 

Do these conclusions hold in cases of other intrinsic emission forms and other values of 
the Lorentz factor discussed above in Sections 4-6? The answer is yes. Presented in Fig. 15 
is the peak energy evolutionary curve arising from a "broad" co-moving pulse of the broken 
power law emission (55) with an unchanged intrinsic spectrum in the decaying phase of the 
pulse, which is seen to possess the most deviation from the marginal decaying curve (59), 
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detected from the peak energy evolutionary curves of all the situations discussed above. We 
find that the conclusions are not violated by all these factors. 

8. Conclusions 

It is known that the power indexes in the relation between the flux and peak energy 
derived from many GRBs are different from that expected by the curvature effect, 7 ~ 3. 
Some possible interpretations for this variation were proposed by authors of [16]: adiabatic 
and radiative cooling processes that extend the decay timescale; a nonuniform jet; and the 
formation of pulses by external shock processes. According to our analysis, the first one is 
unlikely. The cooling time, long or short, is in fact included in the decaying phase of co- 
moving pulses, and the decaying portion of the pulses has little influence on the power index. 
Authors of [15] proposed that the breaking of local spherical symmetry such as a prolate 
or oblate shell geometry would result in a power law index that differs from the spherical 
case. This interpretation is somewhat similar to the second proposal of [16] and might be an 
interesting outlet for the problem. In an investigation on the hardness-intensity correlation 
(HIC) in GRB pulses, authors of [41] found that some pulses exhibit a track jump in their 
HICs, in which the correlation jumps between two power laws with the same index. This was 
interpreted as a signature of the existence of strongly overlapping pulses. This mechanism 
would naturally explain why the index observed in some GRBs could be both larger and 
smaller than what the curvature effect predicts. For example, an upper power law line in 
Fig. 9 jumping to a lower power law line in the figure would give rise to smaller index, while 
a lower line jumping to an upper line would produce a larger index. 

In the relation between the integrated flux and peak energy, the following conclusions 
hold: in the decaying phase of pulses, the two quantities are well related by a power law where 
the index is about 3, being free of the intrinsic emission and the Lorentz factor; the relation 
in the rising phase differs significantly and the overall relation between the two quantities 
varies and shifts enormously, depending on the form or width of the intrinsic emission of 
pulses and on parameters such as the Lorentz factor. 

The pattern of the spectral evolution of the intrinsic emission of a co-moving pulse 
within its rising phase could be observed in the relation between the integrated flux and 
peak energy. For an unchanged intrinsic spectrum, the relation in the rising phase is a 
straight line in parallel with the axis of the flux; for a hard-to-soft intrinsic spectrum, the 
flux decreases with the increasing of the peak energy within this phase; for a soft-to-hard- 
to-soft intrinsic spectrum, the flux generally increases with the increasing of the peak energy 
in the rising phase. 
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Besides the above conclusions, the following could also be reached for pulses arising from 
a relativistically expanding fireball, a) The spectrum of a pulse in its decaying phase differs 
slightly for different intrinsic spectral evolution patterns, and hence it must be dominated 
by the curvature effect, b) An intrinsic soft-to-hard-to-soft spectral evolution within a co- 
moving pulse would give rise to a pulse-like evolutionary curve for the peak energy, c) 
There exists the marginal form of 1/(1 + r/r^) for the peak energy evolutionary curve in the 
decaying phase of pulses, and in many cases the peak energy evolutionary curve well follows 
the form and when the former deviates from the latter it deviates in a reverse-S way. 

Supported by the National Natural Science Foundation of China (No. 10573005 and 
No. 10747001). 
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